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Abstract 

We discuss the importance of flavor ratio measurements in neutrino telescopes, such as 
by measuring the ratio between muon tracks to cascades, for the purpose of extracting 
new physics signals encountered by astrophysical neutrinos during propagation from the 
source to the detector. The detected flavor ratios not only carry the energy information of 
specific new physics scenarios which alter the transition probabilities in distinctive ways, 
but also the energy dependent flavor composition at the source. In the present work, we 
discuss the interplay of these two energy dependent effects and identify which new physics 
scenarios can be distinguished from the detected flavor ratios as a function of astrophysical 
parameters. We use a recently developed self-consistent neutrino production model as our 
toy model to generate energy dependent source flavor ratios and discuss (invisible) neutrino 
decay and quantum decoherence as specific new physics examples. Furthermore, we identify 
potentially interesting classes of sources on the Hillas plot for the purpose of new physics 
searches. We find that sources with substantial magnetic fields 10'^ Gauss < -B < 10® Gauss, 
such as Active Galactic Nuclei (AGN) cores, white dwarfs, or maybe gamma-ray bursts, 
have, in principle, the best discrimination power for the considered new physics scenarios, 
whereas AGN jets, which typically perform as pion beam sources, can only discriminate few 
sub cases in the new physics effects. The optimal parameter region somewhat depends on 
the class of new physics effect considered. 
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1 Introduction 



Even though extraterrestrial high energy neutrino flux has not been detected yet, efforts are 
made to detect neutrinos from far beyond the Sun and the Supernova 1987A. Astrophysical 
neutrinos at high energies open up an entirely new window to infer the properties of pro- 
duction sites as well as neutrinos themselves. The kilometer-scale neutrino telescopes such 
as IceCube [l] and KM3NeT [2] are designed to detect high energy {E > 10^^ eV) neutrinos 
from various astrophysical sources. Among the candidate sources, the most prominent extra- 
galactic ones are Active Galactic Nuclei (AGNs) [3]-[6] and Gamma- Ray Bursts (GRBs) [T], 
see Ref. [s] for a review and Ref. |9| for the general theory governing photohadronic neutrino 
production from astrophysical sources. There are generic theoretical bounds for the diffuse 



neutrino ffux from sources transparent as well as opaque to neutrons 10,11 . The diffuse 



ffux bounds from these sources are being touched by IceCube, see, e.g., Ref. 12 for GRBs 



Apart from that, there is also the possibility to detect point sources for which even the 



optical counterpart is blocked, the so-called "hidden sources" 13 . 

We focus on the possibility to detect the flavor of astrophysical neutrinos. Although neutrino 
telescopes such as IceCube can not detect the flavor of the neutrinos directly, one can, in 
the simplest case, make use of the two kinds of event topologies: tracks (which are induced 
by muons) and cascades/showers (electromagnetic and hadronic, induced by and z/^, 
respectively) to construct an observable track to cascade flavor ratio as a flavor dependent 
quantity. Note that this flavor ratio is deflned in the spirit of the very flrst search for 
cascades from extragalactic neutrino sources, which has been very recently conducted in 



Ref. 14 for IceCube-22. This flavor ratio allows us to infer particle physics properties of 



neutrinos [15^20], extract information on the flavor mixing parameters 21-35 and identify 
such sources |36|j37j. Since the flux normalization drops out of this ratio, it is relatively 
robust with respect to astrophysical uncertainties. We will elaborate on this point further 
in Sec. [HI 

The high energy neutrinos are conventionally expected to originate from the decay of charged 
pions (pion beam source) produced via photohadronic processes (p7) or inelastic (pp) col- 
lisions. From IT ^ fi ^ e decay chain, the flavor composition at the source is given by 
$° : <l>° : <l>° = 1 : 2 : ($° represents sum of neutrino and anti-neutrinos of a given 
flavor) Q However, this picture is over-simplifled and it was pointed out in Refs. [9|[39] that 
energy losses in strong magnetic flelds, which dominantly affect the muons for a pion beam 
source, changes $g : $° : $° to : 1 : at high energies (muon damped source). This 
established the energy dependence of the flavor composition of a given source and one can 
expect a smooth transition from one type of source to another as a function of the neutrino 



energy 40 42], mostly depending on the cooling processes of the intermediate muons, pions 
and kaons (or even charmed mesons, which we do not consider). Other mechanisms, such as 
neutron decay, deflne a new class of sources, and even all special classes may be recovered as 



a function of energy 43 . At low energies, neutron decays may dominate the flux, where the 



^ At neutrino telescopes without charge identification capabiUty, p7 and pp processes are indistinguishable 
since both lead to the same flavor composition 1:2:0 (neutrinos and anti- neutrinos are added). However, 
one may use the Glashow resonance {E ~ 6.3 PeV) initiated by De to discriminate between the primary 



processes which differ in 7r+ — tt symmetry 25 38 
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neutrons are generated in the photodissociation of heavy nuclei or photohadronic interac- 
tions, which leads to 1 : : (neutron beam source) |44] . If the cooled muons pile up at low 
energies, one may even have 1:1:0 from muon decays only (muon beam source) [43]. A 
simple-minded model for implementing the photohadronic interactions in cosmic accelera- 
tors was recently developed by Hiimmer, Maltoni, Winter and Yaguna [HMWY] [43]. In the 
HMWY model, charged pions are produced from photohadronic [pj) interactions between 
protons and the synchrotron photons from co-accelerated electrons (positrons). The pho- 
tohadronic interactions are computed using an efficient state-of-the-art method described 
in Ref. |45|, based on the physics of SOPHIA 46 , and the helicity-dependent muon decays 



from Ref. I41j are included. The toy model relies on relatively few astrophysical parameters, 
the most important ones being the size of the acceleration region (R), the magnetic field 
strength at the source (B), and the injection index (a) which is assumed to be universal for 
protons and electrons/positrons. In the HMWY model, the energy dependent effects men- 
tioned above are taken into account automatically: the synchrotron cooling of all secondary 
species is included, as it is important for the accurate prediction of the fiavor ratios at the 
source. 

The observed neutrino fiavor composition at the detector is in general different from that 
at source due to neutrino flavor oscillations. Standard flavor mixing leads to achromatic 
(energy independent) transition probability for astrophysical neutrinos because of decoher- 
ence of the oscillations. For a pion beam source and standard flavor mixings, one gets an 
interesting prediction of flavor equipartition 1 : 1 : 1 at the detector It should be noted 
that there is no physical basis for this equipartition, it is purely accidental and just a con- 
sequence of the speciflc choice of source and the mixing matrix in the near tri-bimaximal 
form along with the assumption of standard mass-induced oscillations. Hence if we consider 
a source other than pion beam source and/or non-standard propagation effects en-route, 
this prediction of universal ratio changes [l9|. This clearly implies that in order to infer 
new physics one has to be careful, and the mere departure from flavor equipartition can 
not serve as a guaranteed indicator of the speciflc class of new physics in the most general 
situation. The new physics effects studied in the context of high energy astrophysical neu- 

pseudo-Dirac nature of neutrinos 
Lorentz invariance violation 



trinos include neutrino decay 47 51 



of discrete symmetries such as CPT 



51 53 
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violation 
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quan- 



tum decoherence 51,52 54 56 , violation of unitarity of the mixing matrix 57 , coupling of 
neutrinos to dark energy~[58| and non-standard interactions 59 . However in the existing 



studies, the energy dependence of the flavor composition at the source has not been taken 
into account explicitly and many of these studies apply to only a speciflc source type. In 
the present work, we incorporate 



1. Energy dependent effects in the flavor composition at the source by using the HMWY 
model as our toy model, and 

2. Energy dependent new physics scenarios during propagation of neutrinos from source 
to detector. 

We study their impact on the detected flavor ratios for a point source. Most importantly. 



^As long as 



7r/4 and ~ 0, this result is robust and independent of the value of 9i2. 
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unlike the case of standard oscillations which are achromatic, the energy dependent new 
physics effects lead to energy dependent probabilities. Thus energy dependent observable 
flavor ratio is therefore an outcome of an interplay of two energy dependent effects and one 
has to be careful to interpret these results for the new physics searches. In order to elucidate 
this interplay, we explore two representative examples (a) neutrino decay, and (b) quantum 
decoherence. 

The first non-standard possibility considered here is that of neutrino decay over astrophys- 
ical L/E scales, which sheds light on the lifetime of the neutrinos. Among the terrestrial 
experiments, solar neutrino data sets the strongest limit on neutrino lifetime to mass ratio 
/ttl > 10~^ s/eV [47] which implies that decay of high energy astrophysical neutrinos can 
not be ruled out. Astrophysical neutrinos also provide a much higher (by several orders 
of magnitude) neutrino lifetime sensitivity due to the higher L/E ratio involved. Decay of 



astrophysical neutrinos has been studied in Ref. 47 49j. The other new physics scenario 
considered in the present work is that of quantum decoherence effects. Even though to 
date there is no convincing theory of quantum gravity, it is expected to give rise to dis- 
tinctive signatures such as violation of Lorentz invariance, CPT violation and/or quantum 



decoherence 60 at low energies (compared to the Planck scale). Neutrino telescopes are 



particularly well-suited to probe such new physics effects in both atmospheric 61 as well as 



astrophysical 53,54 neutrino fluxes. Note that some of these exotic effects are related due 



to the CPT theorem, for example, CPT violation implies Lorentz invariance violation but 



not vice-versa, and quantum decoherence can give rise to CPT violation 60 . In particular, 
quantum decoherence leads to evolution of pure states to mixed states via interaction with 
the environment of the space-time. The effect of quantum decoherence on astrophysical 
neutrinos has been studied in Refs. I52l[54ll56l . 



The paper is organised as follows. In Sec. [2| we introduce our theoretical flavor ratio 
framework. Then in Sec. |3} we briefly recapitulate the HMWY model and illustrate how 
the flavor ratios of the sources acquire their energy dependence. We introduce our energy 
dependent new physics effects in Sec.|4| where we also discuss their interplay with the energy 
dependent flavor ratios at the source and detector. Finally, we perform a systematic scan of 
regions of the parameter space in Sec. [5| to identify source classes which may be most useful 
for this application. We then summarize in Sec. [6] Note that a more detailed discussion of 
the quantum decoherence model can be found in Appendix |A] 

2 Flavor ratio framework 

We begin by defining the quantities of interest for the present study: the flavor ratio at the 
source, propagation effects, and the flavor ratio at the detector. 
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2.1 Source flavor ratio 



Assuming a negligible amount of Vr at the source the flavor composition is completely 
characterized by the ratio of the electron to muon neutrino flux 

where $e(-^) a^*i "^/K-^) the fluxes of electron and muon neutrinos without propagation 
effects. These correspond to the fluxes directly at the source, apart from an overall normal- 
ization, distance dependence, a possible Lorentz boost, and redshift effects. Since all these 
effects affect the different flavors in the same way, it is convenient to normalize to the fluxes 
at the detector without propagation effects $° The fluxes are typically given in units of 
GeV~^ cm~^ s~^ (point source) or GeV~^ cm~^ s~^ sr~^ (diffuse flux). Note that we always 
refer to the sum of neutrino and antineutrino fluxes in the following, since we assume that 
the detector cannot distinguish between these. 

In the literature, the following source classes are distinguished depending upon the value of 
X: 

Pion beam sources produce neutrinos from charged pion and successive muon decays, 
such as 

7r+ ^ + (2) 
e+ + z/e + z/^ . (3) 

This leads to X ~ 1/2. 

Muon damped sources produce neutrinos from pion decays only, z.e., Eq. (|2]), since the 
muons loose energy by synchrotron radiation efficiently. This leads to X ~ 0. 

Muon beam sources The muons may pile up at lower energies, where muon decays dom- 
inate, see Eq. ([s]). Then we have X ~ 1. Note that the neutrinos from semi-leptonic 
decays of charmed mesons in baryon rich astrophysical environments such as slow jet 
supernovae [ 63] (where pp interactions are prevalent) also lead to high energy neutrinos 
with X ~ 1 1 64]. 

Neutron beam sources The neutrons, produced by photodissociation of heavy nuclei or 
photohadronic interactions, decay into neutrinos by 

n p + e" + z/e , (4) 

leading to X ^ 1. 



•^One may get a Vr component in some exotic scenarios, such as cosmic defects or evaporating black 
holes, typically at energies beyond ^ 10^^ eV. However, the HMWY model does not extend to such high 



energies in most cases because of the Hillas criterion 62 and synchrotron cooling of the protons limiting 
the maximal energy. In addition, one typically expects a diffuse flux from such objects, where as for point 
sources, the angular resolution of the instrument can be used to suppress the backgrounds. Therefore, the 
contribution from these objects can be neglected if point source fluxes are considered. 
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Undefined sources Several processes compete with similar magnitudes, leading to arbi- 
trary X. 

Note that X{E) is, in general, an energy dependent quantity, and these classifications only 
apply to certain energy ranges. A given fiavor composition at source uniquely identifies a 
particular mechanism provided we have some information on the type of source giving rise 
to the neutrino fiux. For instance, a muon beam source can be mimicked by a transparent 
source (via decay of piled up muons at low energies) or a hidden source such as slow jet 
supernovae (via decay of charmed mesons) where different mechanisms are responsible for 
arriving at X ~ 1. In addition, the guaranteed flux of extremely high energy neutrinos 
from Greisen-Zatseptin-Kuzmin (GZK) process (photohadronic interaction of highest energy 
cosmic rays with cosmic microwave background radiation) also yields energy dependent 



source flavor ratios 65,66 . Below 100 PeV, GZK neutrinos fall in neutron beam source 



class while above 100 PeV, GZK neutrinos are more like a pion beam source. 
2.2 Standard propagation effects 

We assume propagation of astrophysical neutrinos in vacuum over distances long enough 
such that oscillation effects become decoherent. On incorporating the propagation effects, 
the flavor flux changes to ^^'^^ for neutrino of flavor zz/j 

^T{E)= J2 Pa,{E)^l{E), (5) 

a=e,fj,,T 

where for standard oscillations, we have 
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Pap{E) = P,^ = V \Upi\' |f/„,|2 , (6) 



i=l 

which is achromatic. Note that Eq. ^ depends only on the mixing angles 9ij contained in 
the mixing matrix elements fZ/jj of the 3x3 Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neu- 
trino mixing matrix which in the commonly adopted Particle Data Group (PDG) parametriza- 
tion is given by 

C12C13 S12C13 sige-^^c^-p' 

U = I -S12C23 - Ci2S23Si3e'^'^P C12C23 - Si2S23Si3e'^'^P S23C13 | , (7) 
S12S23 - Ci2C23Si3e*'^C'P -C12S23 - Si2C23Si3e'^'^P C23C13 

where Sij = sin% and Cij = cos% and Sep is the Dirac type CP phase. Unless noted 



otherwise, we use the following values for the mixing parameters, see, e.g., Ref. 67 : 

sin^ 023 = 0.5, sin^ 612 = 0.318, sin^ ^13 = . (8) 

Thus, we note that flavor equipartition at the detector ^^'^^ : : ^^'^^ = 1 : 1 : 1 is a 

consequence of two inputs: the source ratio X ~ 0.5 (pion beam source) and the use of 
mixing angles being close to the tri-bimaximal form. 
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2.3 Flavor ratio at the detector 



We need to take into account tlie propagation effects, before computing the observable fiavor 
ratios at the detector. Ideally one would want to detect aU the flavors separately, however 
in practice it is not so easy. Several flavor ratios have been constructed in literature to 



distinguish between flavors 20 . The easiest possibility to measure flavor ratios at neutrino 
telescopes requires, apart from muon tracks sensitive to z/^, the identification of cascades [li] . 
These come with a lower statistics and have a higher threshold (about 1 to 10 TeV in 
IceCube). In addition, the neutrino effective area increases much weaker with energy, which 
means that the best statistics may be obtained close to the threshold. However, one expects 
a much better energy resolution for cascades. If we assume that electromagnetic (from z/g) 
and hadronic (from Ur) cascades do not need to be distinguished, a useful observable is the 
ratio of muon tracks to cascades ||23j 

^ _ Pe,iE)XiE) + P,,iE) 

Note that the above formula holds even if unitarity is violated, i.e., Pea + P^J.a + PTa < 1, such 
as for neutrino decay into invisible states. In addition, note that neutral current events will 



also produce cascades, which, in practice, have to be included as background. In Ref. 14 
the most recent IceCube cascade analysis, the contribution of the different flavors for a E~'^ 
extragalactic test flux with equal contributions of all flavors at the Earth was given as: 
electron neutrinos 40%, tau neutrinos 45%, and muon neutrinos 15% (after all cuts). This 
implies that charged current showers dominate and that electron and tau neutrinos are de- 
tected with comparable efficiencies, i.e., that Eq. (|9]) is a good first approximation to discuss 
flavor at a neutrino telescope. The benefit of this flavor ratio is that the normalization of 
the source drops out. In addition, it represents the experimental flavor measurement with 
the simplest possible assumptions. Eq. (joj) implies that the observable flavor flux ratio R 
depends on two energy dependent quantities: X{E) characterizing the energy dependence 
of the flavor composition at the source, and Pai3{E) characterizing the energy dependence 
during the propagation from source to detector. In the absence of energy dependent new 
physics, the probability is independent of energy (see Eq. ^) and the flavor ratio R closely 
follows the energy dependence of the source function X{E). 

Using Eq. ^ in Eq. we have listed the characteristic values of flavor ratios X,R in 
Table [T] assuming standard oscillations and best-fit values of oscillation parameters. In this 
study, we do not consider the uncertainties coming from the oscillation parameters, since 
we expect that the discussed measurements are only feasible with sufficient statistics on a 



timescale when the oscillation parameters are sufficiently limited; see discussion in Ref. 43 
(Fig. 11). 

3 Energy dependent flavor ratios at the source 

The HMWY model describes neutrino production via photohadronic (p7) processes for 
transparent sources (optically thin to neutrons) and includes the cooling of the secondary 
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Source 


X 


R 


Muon damped 





0.64 


Pion beam 


0.5 


0.50 


Muon beam 


1 


0.44 


Neutron beam 


> 1 


0.28 


Undefined 


Other 


Other 



Table 1 : Characteristic values of flavor ratios X and R for the considered source classes for best-fit values 
of oscillation parameters. The ratio R is computed using Eq. and standard oscillation probabilities 
(Eq. (§). 



particles. It can be used to generate neutrino fiuxes as a function of few astrophysical 
parameters. Below we outline the key ingredients of the model and the main results, for 



details see Ref. 43,45 . All of the following quantities refer to the frame where the target 



photon field is isotropic, such as the shock rest frame (SRF). 
3.1 Ingredients of the HMWY model 

The protons and electrons/positrons are injected with spectra oc The maximal energies 
of these spectra are determined by balancing the energy loss and acceleration timescale given 
by 

= ri— , (10) 

with 7] an acceleration efficiency depending on the acceleration mechanism, where we choose 
7] = 0.1 later. If synchrotron losses dominate, the maximal energy is therefore given by 



^-x=y— (11) 

It scales oc m^, which means that the protons are accelerated to much higher energies, and 
oc 1/y/B, which means that strong magnetic fields limit the maximal energies. For each 
particle species, the injection and energy losses/escape are balanced by the steady state 
equation 

QiE) = ^{mN{E)) + ^, (12) 

with tesc{E) the characteristic escape time, b{E) = -Et^^]^ with t^J,,{E) = -1/EdE/dt the 
rate characterizing energy losses, Q{E) the particle injection rate [(GeVscm^) ] and N{E) 
the steady particle spectrum [(GeVcm^) ^]. For all charged particles, synchrotron energy 
losses and adiabatic cooling are taken into account. In addition, unstable secondaries, i.e., 
pions, muons, and kaons, may escape via decay. As a consequence, for pions, muons, and 
kaons, neglecting the adiabatic cooling, the (steady state) spectrum is loss-steepend above 
the energy 
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Qe(E) [->(N;(E)y^ Q,(E) [->(N,(E)^^N7ijy-] Q^{e) 



photohadronics 
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N.(E) 



weak decays 
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from fj* 



from // from /i 



Figure 1: Flowchart describing the neutrino production in the HMWY model. The functions Q{E) 
denote (injection) spectra per time frame [(GeVscm^) ^] and N{E) steady spectra [(GeVcm^) ^] derived 
from the balance between injection and losses or escape. Dashed arrows stand for solving the steady state 



differential equation Eq. (12), the horizontal line "photohadronics" to solving Eq. (16 1 for all interaction 



types. Figure taken from Ref. 43 . 



where synchrotron coohng and decay rates are equaL One can read off this formula that the 
different secondaries, which have different masses m and rest frame lifetimes tq, will exhibit 
different break energies Ec oc ^/vtv' jr^ which solely depend on particle physics properties 
and the value of B. 

While being accelerated, the electrons loose energy into synchrotron photons, which serve 
as the target photon field. Charged meson production then occurs via 

p + 7-)-7r+p', (14) 
p + 7^i^+ + A/S, (15) 

with these synchrotron photons, where the leading kaon production mode is included and p' 
is a proton or neutron. In addition, two- and multi-pion production processes are included 



not listed here), see Ref. 45 for details. The injection of the charged mesons is computed 



from the steady state proton Np{Ep) and photon N^{e) spectra with 45 

oo oo 

Q,{E,) = j ^ Np{Ep) j deN,{e)Rb{x,y), (16) 

Eb ^ Hh'^P 

2Ep 

with X = Eh/Ep the fraction of energy going into the secondary, y = {Epe)/mp (directly 
related to the center of mass energy), a "response function" Ri,{x,y) (see Ref. |45j), and eth 
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the threshold for the photohadronic interactions (in terms of photon energy in the proton 
rest frame). The weak decays of the secondary mesons, such as Eqs. (|2| and ([s]), are 
described in Ref. [4l]. These will finally lead to neutrino fluxes from pion, muon, kaon, 
and neutron decays. The whole procedure is illustrated in Fig. [T} the finally obtained eight 
neutrino fluxes are shown as shaded boxes. In the following, we will sum over the fluxes 
from all polarities. For instance, if we refer to "z/^ from muon decays", we mean the sum 
over the from fi~ and the z/^ from /i"*" in Fig. [l] 



One can clearly see from Eq. (16), that only the product of the proton and photon (and 
therefore electron) density normalizations enters the final result, as long as cooling processes 
implicitly depending on proton-electron ratio, such as inverse Compton scattering or photo- 
hadronic cooling, are sub-dominant. This is, of course, a key simplifying assumption which, 
in general, limits the applicability of this model. However, it leads to the simplest possible 
model which includes the effects relevant for the flavor ratios, whereas additional cooling 
processes affecting the protons and electrons only change the shape of the neutrino spectra 
(including maximal energies), but not the flavor ratios as a function of energy directly. In 



the flavor ratios, the overall normalization from Eq. (16) cancels. In that sense, the flavor 
composition might be the most robust prediction one can make for a source, since it is very 
little affected by astrophysical uncertainties as long as the dominating cooling process of 
the secondaries (pions, muons, kaons) is synchrotron emission, and the dominant escape 
process is decay. These processes depend, however, on the magnetic field in the source B, 
a quantity, which is not directly accessible. 

In the following, for the sake of simplicity, we assume that we can estimate the parameters 
of the source by utilizing the multi-messenger connection (such as gamma-ray observations). 
For example, one can estimate the magnetic field B from energy equipartition, the injection 
index from the spectral shape, and the size of the acceleration region R from the variability 
timescale. Then one can predict the flavor ratios as a function of these parameters. Con- 
versely, flavor ratio measurements may provide a direct handle on astrophysical quantities 



such as the magnetic field 43 . We use a = 2, since we do not obtain strongly qualitatively 



different results in the range 2 < a < 3. 

3.2 Hillas criterion and source classification on the Hillas plot 

The parameters R and B can be directly related to the Hillas plot, see Fig. |2j left panel, 
for an example. In order to confine a particle in a magnetic field at the source, the Larmor 
radius has to be smaller than the extension of the acceleration region R. This can be 



translated into the Hillas condition for the maximal attainable energy 62 



£^max [GeV] ~ 0.03 -ri-Z-R [km] ■ [G] , (17) 

where Z is the charge (number of unit charges) of the accelerated particle, B is the magnetic 
field in Gauss, and rj can be interpreted as an efficiency factor or linked to the characteristic 
velocity of the scattering centers. One complication in this type of figure is that R and B in 



Eq. (17) are potentially given in the SRF, whereas -Emax is to be determined in the observer's 
frame. If the particles are accelerated in a relativistically moving environment, such as in 
a GRB fireball, this assignment is not trivial anymore, and the relatively large Lorentz 
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Figure 2: Left panel: Possible acceleration sites in Hillas plot as a function of R and B (version adopted 
from M. Boratav). Right panel: Classification of sources for injection index a = 2 in this plot (see main 
text). Some points from left plot are shown for orientation, as well as two new points (12 and 13) are defined 
for later reference. Figure taken from Rcf. ,43j . 

boost r of the acceleration region into the observer's frame must be taken into account. 
We stick to the interpretation of R and B in the SRF, which means that (for 77 = 1) the 



condition in Eq. (17) depends on the Lorentz boost of the source. This is illustrated by 
showing Eq. (17) for several selected Lorentz boosts in Fig. [2| left panel, for protons, ?7 = 1, 
-f'max = 10^° eV. In the following, we will use F = 1. However, it should be noted that if the 
source is significantly boosted, the neutrino energies will have to be increased by F, and the 
(cascade) threshold of the neutrino telescope may be passed. We neglect redshift effects on 
the other hand. 

In Fig. [2| right panel, the main result of Ref. |43:| is shown for a = 2: Here the sources 
are classified into the categories from Sec. [2] as a function of R and B. Note that a given 
source may fall into different categories in different energy ranges. For example, the "muon 
beam — )■ muon damped" region means that a muon beam is found for low energies, whereas 
a muon damped source is found for higher energies. The classical "pion beam — )■ muon 
damped" category is found for large enough magnetic fields and large enough R as well, 
where the transition between adiabatic and synchrotron cooling leads to an additional break 
in the proton spectrum. 

We illustrate the classification of the sources for test point (TP) 13 from Fig. [2] (right panel), 
see Fig. [3} where the neutrino fluxes and flavor ratio at the source are explicitly shown for 
an injection index, a = 2. One can easily see in the left panel that the different cooling 
and decay timescales of the secondaries lead to a hierarchy of the breaks, as described in 



Eq. (13), an effect which is similarly present in GRBs (see, e.g., Refs. 41,68]). The Uf, flux 
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10"- 10° lO" lO"* 10* 10* 10"^ 10° 10- lO** 10*^ 10* 10"^ 10° 10^ lO'' lO" 10* 

Energy [GeV] Energy [GeV] Energy [GeV] 



Figure 3: Muon (left panel) and electron (middle panel) neutrino fluxes at the source, and X (right 
panel) as a function of energy for test point (TP) 13 from Fig. [2j right panel. The contributions of the 
individual neutrino fluxes from different parents are shown as well. The gray-shaded regions mark the 
energy window used for the source classiflcation (see main text). Figure adopted from Ref. [43j . 

mostly comes from muon decays (middle panel), which means that it will roughly follow 
the muon break as in the left panel, whereas from pion and kaon decay (left panel) 
extend to higher energies. As a consequence, the flavor ratio X changes from 0.5 to 0.0 
in the right panel at about 10^ GeV. The neutrinos from neutrons, which are produced in 
the photohadronic interactions, lead to a z/g flux dominating at low energies. However, in 
that range the total flux is already too low compared to the peak and also drowned in the 
atmospheric neutrino background. 

Of course the above identification of different source types applies only to an energy range 
where the flux of at least one of the flavors is large enough, to give reasonable statistics. 
Since the region close to peak in E^(^^^ will contribute mostly to the event rates, at least in 
energy ranges where the neutrino effective area is proportional to E"^, we define an energy 
window which captures the "upper two orders of magnitude" in the flux. That is, we 
compute (for each point) the maximal flux in and derive the energy range where the 

flux of any flavor is at least 1% of the respective maximal flux (in ii^^$^). Only this energy 
range is considered, since otherwise the flavor ratio may be ill-defined. This energy range 
is marked as shaded regions in Fig. |3j In the right panel, one can also see that the flavor 
ratio is ill-deflned above 10^ GeV, where the fluxes are negligible. In addition, note that it 
is required that a speciflc flavor ratio category is found over at least one order of magnitude 
in energy within the chosen energy window. 

4 Energy dependent new physics scenarios 

In presence of new physics, the flavor ratio R not only depends on energy dependent flavor 
composition of the source X{E), but also on the new physics induced energy dependent 
terms in the probability, see Eq. ^ . We consider two speciflc examples of new physics here 
which can give rise of energy dependent effects in opposite energy regimes. Whereas neu- 
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trino decays are mostly present at lower energies (the lifetime is Lorentz-boosted), quantum 
decoherence may plausibly cause the strongest effects at higlier energies. The energy de- 
pendence of source is contained in X{E) as described in Sec. [sj while the energy dependent 
new physics effects lead modification of probabilities relevant for astrophysical neutrinos as 
described below. 



4.1 Neutrino decay 



Neutrino decay is usually described by an energy dependent damping of the overall oscilla- 
tion probability 47 and Eq. ^ gets modified to 



Pa/3 = XI|t^^i|'|?7«^r A(^) with D,{E) = exp (-a~] , (18) 
i=i ^ ^ 

as the damping coefficient^ Here dj = mi/r^ with is the rest frame lifetime for mass 
eigenstate z/j. Typically the neutrino lifetime is quoted as rf /uii since rrii is unknown. From 
the exponential factor, the neutrinos decay if 

^=«r^ < lO^^^seV-^ (19) 
rrii Mpc h 

Thus, neutrino telescopes can probe lifetimes of the order of ~ 10^ seV~^ for L ~ 1 Mpc 
(~ 10^^ s) and E' ~ 1 TeV. In what follows, we use the same lifetime for all mass eigenstates, 
rf /nii = T^/m for the sake of simplicity. 

In view of the rather weak direct limits on neutrino lifetime and the proposed unparticle 
scenarios, one can have many different mechanisms of neutrino decay. In what follows, we 
do not assume a particular decay scenario but consider the general case j49]. For neutrino 
decay, one typically distinguishes: 

Visible decays One of the decay products is visible (to the neutrino detector), typically 
a lighter active neutrino mass eigenstate. 

Invisible decays The decay products are invisible (to the neutrino detector), such as ster- 
ile neutrinos or unparticles. 



In addition, we have 



Incomplete decays The mass eigenstates have decayed partially, i.e., < Di < 1. 
Complete decays The mass eigenstates have decayed completely, i.e., Di — 0. 



In our framework Eq. (18), we can describe complete or incomplete invisible decays. For a 
framework for visible complete decays, see Ref. 49 . For the most general framework, see 
',7l]. Note that if only a single mass eigenstate is stable in complete decays, R does 
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Figure 4: Flavor ratio i? as a function of X for all complete decay scenarios. Black disks refer to stable 
mass eigenstates, white disks to unstable mass eigenstates. Different sources classes as a function of X are 
marked. 



not depend on X{E) since probabilities factorize in the source-dependent and detector- 



dependent parts; c/., Refs. 15,49. In general, any decay scenario is characterized by 
stability properties of active states which gives in all 2^ = 8 scenarios for invisible decays, 
since each mass eigenstate can be either stable or not. We do not consider the case of only 
unstable states, since no signal will be observable then. 

We show in Fig. [4] the flavor ratio i? as a function of X for all complete decay scenarios. 
Black disks refer to stable mass eigenstates, and white disks to unstable mass eigenstates. 
One can read off this figure, that for the pion beam, four scenarios can be clearly separated, 
whereas three scenarios are almost degenerate. If, however, the information from a muon 
damped sources is added, such as at high energies, exactly these three scenarios, in principle, 
separate. The muon beam has a similar effect in the opposite direction of X. Therefore, it 
is crucial to discuss sources which contain information from different parts of X. 

We show in Fig. [s] two anomalous sources described in Sec. [sj TP 3 (upper row) and TP 13 



(lower row) from Fig. |2} In the left panels, the energy dependent damping coefficients Da 



(for aL = 10^ GeV 
at the source X{E) 



and Db (for aL = 10 GeV) are shown, as well as the flavor ratio 
^ One can read off these panels, that decays become effective below 



'^In the spirit of Ref. [69], see Eqs. (2.3) and (2.4), which contain different damping effects. For astro- 
physical neutrinos, we can only probe Da = Di. 

^Note that the chosen decay rates correspond to lifetimes of about 0.1 seV~^ (Da) and lO^"^ scV^^ (Db) 
for L = 1 Mpc using /m ~ . These lifetimes are at least one order of magnitude above the current 
direct bounds (depending on the mass eigenstate) except from the SN 1987A bound (see, e.g., Ref. [20] and 
references therein). It is, however, strictly speaking not trivial to assign the SN 1987A bound to Ui or ^2, 
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Figure 5: Left panels: Energy dependence of the damping coefficients for aL — 10^ GeV {Da) and for 
aL = 10^ GeV (Db), as well as for the function X{E) (dashed curves). The upper row corresponds to TP 3 
from Fig. [2] (AGN nuclei), the lower row to TP 13. Middle and right panels: Energy dependence (solid 
curves) of the flavor ratio R for seven different decay scenarios. Black disks refer to stable mass eigenstates, 
white disks to unstable mass eigenstates. The middle panels are plotted for aL = 10^ GeV {Da) and the 
right panels for aL = 10* GeV {Db), where the same decay rate is chosen for all unstable mass eigenstates. 
The region where the flux is relatively large is shown as shaded region. The dotted curves refer to complete 
decays. 



10^ GeV (Da) or 10*^ GeV (Db) for our choice of parameters. Well above these energies, 
the mass eigenstates are practically stable (the lifetime increases with the Lorentz boost of 
the particles). The curves of X{E) clearly indicate that both the sources TP 3 and TP 13 
(upper left and lower left panel) perform as muon-damped sources at high energies for which 
X 0. At lower energies, TP 3 is undefined, whereas TP 13 performs, within reasonable 
uncertainty, as pion beam. 

In the middle and right panels of Fig. [sj the flavor ratio R is shown for seven different decay 
scenarios (solid curves). The middle panel corresponds to while the right column to 

given the low statistics; the problem is that a flavor eigenstate {De) is measured which contains significant 
contributions from ui and V2- The only conclusion from that observation is that it is unlikely that both vi 
and V2 are unstable at that energy-distance scales, which means that the upper curves in Fig. [5] (middle 
and right panels) are basically excluded. 
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Db- As mentioned before, the same decay rate is chosen for all unstable mass eigenstates. 
Note that X{E) and R{E) have opposite behavior as a function of E. Let us first focus on 
plots corresponding to Da (middle panel). From the left panel for TP 3 and TP 13, we can 
read off that the Da and X curves intersect at about the same energy, which means that 
neutrino decays are only effective in the undefined (TP 3) or pion beam (TP 13) energy 
range below about 10^ GeV. From Fig. |4| we know that this means that some of the decay 
scenarios cannot be clearly separated, which we indeed find in the lower energy range of 
the middle columns of Fig. [sj For E ^ 10^ GeV, however, all mass eigenstates are stable, 
which can be seen from the deviation of the solid curves with respect to the dotted ones 
which refer to complete decays (independent of energy). If all mass eigenstates are stable, 
the same physics scenario is recovered in all cases, and no new information is added. This 
is different for the higher decay rate in the right column: here the muon damped nature of 
the source allows for a separation of the different curves in energy range E ~ 10^ — 10* GeV, 
where both the test points evolve into muon-damped source type, but the mass eigenstates 
are not yet stable. In this range, the degenerate curves separate out. Finally, Fig. |5]is an 
example for the interplay between the source characteristics X{E) and the decay-type new 
physics scenario. 



4.2 Quantum decoherence 



The impact of quantum decoherence on the flux of high energy astrophysical neutrinos has 
been studied in Refs. 51 ,52,54,56 . We show in Appendix [A| how and under what conditions 



the oscillation probabilities presented for three neutrino flavors in this section are obtained. 
Under the assumption that oscillations average out over astrophysical distances, we can 
express the transition probability as a function of only two non-zero decoherence parameters 
^ and 6, 

Po., = l + \iUl, - Ul,){Ul - U%)D^ + \{Ul, + Ul, - 2Ul,){Ul, + U% - 2Ul,)Ds , (20) 

where D^f and Ds are the damping factors (corresponding to the eigenvalues of A3 and As, 
respectively, of the decoherence matrix described in Appendix |A]) given by 

D^{E) = exp (-2 kLE'^) . (21) 

Here Df^[E) parameterizes effects due to quantum decoherence and Uai are the elements 
of the standard neutrino mixing matrix. Here n carries the energy dependent imprint of 
a specific model. In the literature, n = —1,0,2 have been used (see also Ref. [69j). In 
principle, ^ and 6 can take different values, however, in what follows we will assume the 
same energy dependence for the two parameters. 



From Eq. (20), we note that quantum decoherence scenario always leads to flavor equipar- 
tition (1 : 1 : 1) as we take the asymptotic limit L — > 00 if both 5 > 0, whereas in the 
limit 5 — 0, we recover the standard oscillation result (Eq. (16))). Physically, the meaning 



of flavor equipartition (1 : 1 : 1) is that neutrino flavor gets completely randomized due to 

^Note that, k is used to denote 4', S and compared to Appendix|A) we incorporate the energy dependence 
exphcitly here, i.e., we replace "if ^i?", S — 6E^. 
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interaction of system with the environment and all the flavors get equally populated. This 
result does not depend on the particular type of source or choice of mixing parameters. 
However, as noted before flavor equipartition was also obtained in the standard oscillations 
(without decoherence) specifically for the pion beam source for the best-fit values of mixing 
angles. This coincidence is purely accidental and there is no physical reason for this coinci- 
dence. In order to infer effects due to quantum decoherence therefore the pion beam source 
is rendered useless. However, any of the anomalous sources (which do not lead to acciden- 
tal flavor equipartition) are thought to be good candidates for the study of decoherence 
vis-a-vis standard oscillations. For instance, it is shown that detection of galactic electron 
anti-neutrino beams at IceCube will lead to a major improvement in sensitivity to quantum 



decoherence effects 52 , 54 , 56 . We find such a source at low energies or the extreme high-B 
region in Fig. |2} However, the neutrino energies are typically too low in these cases to have 
flavor identification in the neutrino telescopes (the extremely high magnetic field leads to 
proton synchrotron losses such that it is difficult to accelerate to high energies). 

For a fixed value of n (which characterizes a particular model) , we choose the following four 
cases: 

1. (5 = \E' = (standard oscillations) 

2. 5 = $ > 

3. 5 = and ^ > 

4. ^ = and 5 > 

Since n = leads to an energy-independent effect, we do not consider this case. Therefore, 
we study the following two cases: 

(a) n = 2: This corresponds to a string inspired model [60], 



D^{E) = exp {-2kLE'^) . (22) 

Here the quantum decoherence effects are present at high energies. From the expo- 
nential factor, setting 2kLE'^ ~ 0{1) we have decoherence for 

< 2 X 10^^ ( GeV . (23) 

Mpc VTeV J ^ ' 

This implies that astrophysical neutrinos with E = 1 TeV and L = 1 Mpc allow us to 
probe K, ~ 10~^^ GeV~^ for this case. 



(b) n = — 1: This corresponds to a Lorentz invariant model 60 



D^{E) = exp ( -2k^ ] . (24) 
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Figure 6: Flavor ratio i? as a function of X for all complete decoherence scenarios {L large). Black 
disks refer to coherent parameters {k — 0), white disks to decoherent parameters (k > 0). Different sources 
classes as a function of X are marked. Left figure is for our best-fit values of the parameters (sin^ 623 = 0.5), 
while the right figure is for sin^ 6*23 0.4. 



Here the quantum decoherence effects are present at low energies. The spectral sig- 



nature is similar to decay (Eq. (18)). From the exponential factor, setting 2kL/E 
C(l), we have decoherence for 



< 2x 10^5-^^ GeV-2, 
Mpc E 



(25) 



This implies that astrophysical neutrinos with E = 1 TeV and L = 1 Mpc allow us to 
probe K ~ 10~^^ GeV^ for this case. We note that even though decay and decoherence 
are two distinct physical processes, non-zero value of decoherence paramters \E', 6 can 
be used to set a limit on the neutrino lifetime, from Eq. (25) we get = r^/m ~ 
10^ seV for E = 1 TeV and L = 1 Mpc. 

We have introduced four scenarios for decoherence, since any decoherence parameter (\E' or 
5) can be switched on or not. We show in Fig. [g] the flavor ratio i? as a function of X 
for all complete decoherence scenarios, i.e., the limit of large enough L. Black disks refer 
to coherent parameters (k = 0), white disks to decoherent parameters {k > 0). As for 
the decay, we have chosen 5 = \E' in the case where both parameters are decoherent. First 
of all, one can read off from this figure (left panel for maximal atmospheric mixing) that 
(5 = \Ef > leads to R = 0.5, irrespective of the value of X, since universal flavor mix is 



a generic prediction of the decoherence scenario. This can be seen from Eq. (20). Even 



more surprising is the fact that all four scenarios have the same value of R as for complete 
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Figure 7: Left panels: Energy dependence of the decoherence coefficients for case (a) {Dk,{E) — 
exp {-2nLE^)) for nL = IQ-^^ GeV"^ (Dc) and for nL = IQ-^^ GeV"^ (Do), as well as for tlie function 
X{E) (dashed curves). The upper row corresponds to TP 3 from Fig. [i] (AGN nuclei), the lower row 
to TP 13. Middle and right panels: Energy dependence (solid curves) of the flavor ratio R for the four 
different decoherence scenarios in case (a). Black disks refer to coherent parameters (k = 0), white disks to 
decoherent parameters (k > 0). The middle panels are plotted for kL — 10^^^ GeV^^ {Pc) and the right 
panels for kL = W)-^^ GeV"^ {Do)- The region where the flux is relatively large is shown as shaded region. 
The dotted curves refer to complete decoherence. 



decoherence case (case 2) for the pion beam. This is actually due to the conspiracy between 
the use of best-fit values and the value of X which leads to vanishing of terms dependent 
on the mixing matrix elements in Eq. (20) for other three cases, including that of standard 
oscillations (see Appendix ffl). If we change one of the mixing angles (sin^^23), the curves 



indeed separate out for the pion beam source (see Fig..[6| right panel). The biggest splitting 
is obtained for the neutron beam, whereas the scenarios can also be distinguished for the 
muon beam and muon damped source, in principle. In these cases, the curve \Ef = 5 > 
deviates strongest from the standard case = 5 = 0. Note that compared to neutrino 
decay, all the different scenarios are relatively close to each other, which means that in 
practice it may be extremely difficult to disentangle them. 

We show in Fig. [7| the four possible scenarios for case (a) (n = 2), which are qualitatively 
different from decay because here decoherence enters at high energies. This figure is sim- 
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ilar to Fig. |5| with the white (black) disks referring to decoherence (no decoherence) from 
the specific parameter. The middle and right panels correspond to two different values of 
the decoherence parameters Q As expected from the energy dependence of decoherence 
induced damping, the main effects are present for high energies, where the scenarios dif- 
ferentiate thanks to the muon-damped behavior of both the test points. In these cases, 
they approach the "full decoherence" curves (dotted). Comparing Dc (middle panels) and 
Dd (right panels), Dc catches a part of the transition region between the different source 
types. However, since the scenarios are quite alike there, this does not help. The decoher- 
ence scenario Dd, however, leads to effects practically beyond the peak of the spectrum. 
Therefore, kL = 10~^^ GeV~^ is about the maximal sensitivity one can achieve. Comparing 
TP 3 (upper row) with TP 13 (lower row), one can see from the dotted curves that one 
could, in principle, also use the range E <^ 10^ GeV in the upper middle panel for the case 
of complete decoherence. However, the decoherence effects do not show up at such low 
energies for the chosen parameter values. 

In Fig.jsj the case (b) {n = —1) is shown for comparison. In this case, the energy dependence 
is similar to decay, and the effect of decoherence is present at lower energies. As for decay, 
it is very important here that the decoherence effect is large enough to catch a part of the 
muon-damped range. However, compared to decay, none of the scenarios are distinguishable 
at low energies, and one has to rely on the high energy information. In case (b) (compared 
to case (a)), the range for the decoherence parameters for which the scenarios (for De and 
Dp) can be distinguished, is small. For example, in the middle panels [Dp) the decoherence 
parameters are not large enough to produce observable effects in the muon-damped range. In 
the right panels [Dp), there is some energy range where the scenarios can be distinguished. 
At high energies, the muon-damped fixed point is then approached again by all scenarios, 
which is similar to decay. 

5 Potential sources on Hillas plot for new physics searches 

Here we discuss which source classes may be especially useful for the purpose of extraction 
of energy dependent new physics scenarios. We do not take into account instrumental or 
mixing parameter uncertainties. We also do not perform a statistical treatment, but we 
comment on statistics at the end of this section. We study where to look for possible new 
physics effects from the theoretical point of view, i.e., under ideal set of assumptions. 

As method, we count the number of our (seven) decay or (four) decoherence scenarios which 
we can uniquely determine for any energy around the spectral peak at a specific tolerance. 
That is, for each pair of parameters {R, B}, it is tested if the energy dependent flavor ratio 
can be distinguished from all other scenarios at least at some energy value (which can be 

^For a distance of i = 1 Mpc, the chosen values for case (a) (n — 2) correspond to k ^ 10~^° GeV~^ 
for Dc (kL = 10-12 GeV"2) and k ~ 10"^'' GeV-^ for Dd (kL = lO^^^ GeV-^) respectively. The chosen 
values for case (b) (n = -1) correspond to /c ~ lO'^^ QeV^ for De [kL = 10^ GeV) and k - lO^^o Q^y'^ 
for Dp {kL — 10® GeV) respectively. Specifically, for case (a), that is about 10^'* times more sensitive 
compared to terrestrial long baseline experiments, while for case (b) the sensitivity is about 10^ times 
higher (see Ref. (60]). The reason for higher sensitivities is the much longer distances, and, in case (a), also 
the higher energies involved. 
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Figure 8: Left panels: Energy dependence of the decoherence coefRcients for case (b) {D^{E) — 
exp (-2k|)) for nL = 10^ GeV [De) and for kL = 10^ GcV [Dp), as well as for the function X{E) (dashed 
curves). The upper row corresponds to TP 3 from Fig. [2](AGN nuclei), the lower row to TP 13. Middle 
and right panels: Energy dependence (solid curves) of the flavor ratio R for the four different decoherence 
scenarios in case (b). Black disks refer to coherent parameters (k = 0), white disks to decoherent parameters 
(k > 0). The middle panels are plotted for kL = 10^ GeV {De) and the right panels for kL — 10* GeV 
(Dp)- The region where the flux is relatively large is shown as shaded region. The dotted curves refer to 
complete decoherence. 

different for different scenarios). Then the number of such unique scenarios is counted. 
Note that compared to Fig. [2| where the source is to be found over at least one order of 
magnitude in energy, we also allow for smaller energy windows. In addition, note that this 
approach is conservative from the theoretical point of view, because even though a scenario 
may not be uniquely established, some other scenarios may be excluded. For the numerical 
analysis, we need to impose some numerical uncertainty A on R (tolerance). For example, 
for A = 10% and R = 0.5 for the encountered (true) scenario, only scenarios with R < 0.45 
or i? > 0.55 can be excluded. 

In Fig. |9| we show the number of decay scenarios (color coding: see legend) which can be 
uniquely determined as a function of R and B for a = 2. In the left panel, the tolerance 
A = 10%, in the right panel, A = 30%. The lighter the color, the more scenarios can, 
in principle, be uniquely determined. The (light) spikes at the bottom {R ~ 10^^ km. 
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Figure 9: Number of decay scenarios (color coding: see legend) which can be uniquely determined for 
the model from Ref. [43J on the Hillas plot. In the left panel, the tolerance A = 10%, in the right panel, 
A = 30%. Here aL = 10^ GeV (Db) and a = 2. See main text for details. 



B ~ 10^^ Gauss) come from a contribution from neutron decays. For tlie (dark) spike on 
the right sides [R ~ 10^^ km, B ~ 10^^ Gauss), the energies around the spectral peak are 
already too high, i. e., the neutrinos from these sources are practically stable. If the tolerance 
is low enough, obviously all scenarios can be identified uniquely if the magnetic field is in 
the range to allow for both the pion beam and the muon damped part at high energies; c/.. 
Fig. [5j If only one type of source is dominant, such as for the pion beam region in Fig. |2l not 
all scenarios can be disentangled. For the tolerance A = 30%, of course, fewer new physics 
scenarios can be identified. However, it is clear from Fig. [9} that the potentially interesting 
sources have substantial magnetic fields 10^ Gauss ^ B < 10^^ Gauss, which are, however, 
not too high to lead to the entire domination of the neutron beam or muon damped part. 
For example, for our TP 3 (AGN cores), depending on the tolerance, three to six scenarios 
can be uniquely determined, which is already very good. Other good candidates might 
be white dwarfs and gamma-ray bursts (not shown here, since they are not described well 
within our toy model). 

We show similar analysis for decoherence scenarios (a) and (b) in the left and right panels 



of Fig. [10} respectively. Here we have chosen a much smaller values for the tolerance 
(A = 5%), since all scenarios are very close to each other. Since the decoherence scenarios 
are identical for the pion beam source and split up in equally separated values of R for the 
muon damped source, there is practically no qualitative effect of the number of scenarios: 
either all or none can be identified. For both decoherence scenarios (a) and (b), the pion 
beam source alone is not sufficient, as we discussed earlier, which means that substantial 
magnetic fields are required. As a qualitative difference between scenarios (a) and (b), (a) 
requires high energies for the effect due to decoherence to be present, whereas the signature 
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Figure 10: Number of decoherence scenarios (color coding: see legend) which can be uniquely determined 
for the model from Ref. [43] on the Hillas plot. In the left panel, decoherence scenario (a) is shown for 
kL = 10~^^ GeV~^ (Dc), in the right panel decoherence scenario (b) for kL = 10® GeV {Dp). Here a = 2 
and a tolerance on the flavor ratio of 5% is assumed. See main text for details. 

of (b) is strongest at low energies. Since proton synchrotron losses in the magnetic field 
limit the maximal energy in the considered model, the magnetic field cannot be too strong 
in case (a), whereas the optimal region for (b) is basically consistent with the non-pion 
beam region in Fig. [2j Note, however, that for too large B, the neutrino energies may not 
even be above the cascade threshold in an experiment such as IceCube (unless the source 
energies are highly Lorentz-boosted) . In addition, in case (a), magnetic fields as low as 
10 Gauss may be useful if i? > 10^^ km, because then high neutrino energies are obtained 
while magnetic field effects can be still observed. In both cases, TP 2 and TP 3 lie in the 
optimal region, similar to the decay example above. 

In principle, one can also look into the parameter ranges in the individual scenarios which can 
be constrained. We find that decay and decoherence (b) prefer high magnetic fields, where 
the muon damped (and neutron beam) ranges are found at relatively low energies. However, 
again, note that in this case the event rates are surely extremely low, and that the energies 
may not even be above the cascade threshold in an experiment such as IceCube. On the 
other hand, decoherence (a) scales basically with the maximal energy allowed in the model, 
i.e., lower B and larger R are preferred. However, as discussed, B cannot be too low (lower 
than about 1 Gauss), since the effect cannot be seen for the pure pion beam. Considering 
the optimal tradeoff between high enough neutrino energies and large enough new physics 
effects, probably the lightest region in Fig. |9} left, represents the parameter space of the 
most useful sources for new physics effects present at low energies, and the light region in 



Fig. 10 (left), for effects present at high energies. In the range 10^ Gauss ^ B < 10^ Gauss, 



most distinct scenarios may be accessible. 
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The discussion in this section has been performed without dedicated statistical study, since 
it is not clear which sources will provide signals and how large these could be. The number 
of events N can be estimated by the convolution of the flux with the exposure (neutrino 



effective area times time), see, e.g., Ref. 72 for the effective area of IceCube-40 for dif- 



ferent source declinations. In the absence of backgrounds, which may be used as a first 
approximation for point sources using the angular resolution for background suppression, 
the flux limit scales inversely with A^. At the 90% CL, 2.4 events are still compatible with 
the non-observation of a flux in about one year of IceCube-40 data. Taking into account 
that the effective area of IceCube-86 is about a factor of three to four higher at 100 TeV 



(see Ref. 73 , Fig. 5) and assuming a total exposure of ten years, about 100 (muon track) 
events in total can still be expected in the full-scale experiment if the current bounds are 
saturated. The effective area for cascades, which are dominated by z/g and Uj. is, however. 



about a factor of ten lower at 100 TeV 14 , 73 , where it somewhat helps for statistics that we 



use the fluxes from both and Ur and do not require an unfolding of the z/g events, as, e.g. 



in Ref. 41 . For flavor equipartition at the detector, this means that for 100 muon tracks 
only 10 cascades may be expected, and that cascades may constrain the statistics. This 
however does not apply to all cases. For instance, if only mi is stable (see also Ref. |41|), 
R ~ 0.2 implies that considerably more cascades than muon tracks are expected, up to 50 
cascades are compatible with the current (muon track) bounds, corresponding to a statis- 
tical error of R of only 14%. Of course, a unique scenario identification is very difficult. 



but some scenarios may be excluded then (see, e.g., Ref. 49 for a more detailed statistical 
discussion). In addition, for detectors such as the DeepCore array, more muon tracks are 
fully contained, which means that this discrepancy should be smaller. In Ref. 41 , it was 
also pointed out that the detector response for muon tracks and showers depends on the 
spectral shape, since the energy dependence of these two event classes is different. It may 
be regarded as a strength of our approach that we can not only predict the flavor ratio, but 
also the spectral shape. Finally, mixing parameter uncertainties have not been taken into 
account. However, we expect that these will be reduced by future measurements on the 
time-scale discussed here to a level at which they are only a sub-dominant contribution to 



the total uncertainty, see Fig. 11 in Ref. 43 
6 Summary and discussion 

We have discussed the role of flavor ratio measurements at neutrino telescopes to decipher 
potential new physics effects during the propagation of neutrinos. While in the literature, 
flavor ratios at the source are typically considered as constant numbers for particular source 
classes (pion beam, muon damped, muon beam, neutron beam), we have incorporated energy 
dependence in these flavor ratios at the source. We assume that these flavor ratios can be 
predicted well if the astrophysical parameters of the source are known, such as size of the 
acceleration region R, magnetic fleld B of the source, and injection index a. In many cases it 
is possible to estimate these parameters from the observation of the gamma-ray counterpart, 
such as by time variability, energy equipartition, and the spectral shape, respectively. 



We have used the HMWY toy model 43 for the flux computation, which can describe rel- 



atively wide parameter regions using as few assumptions as possible. In this case, charged 
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mesons are produced from photohadronic interactions between protons and the synchrotron 
photons of co-accelerated electrons/positrons. We have emphasized that the flavor ratios 
for a specific energy do not depend on as many assumptions as the spectral shape. Basi- 
cally, they are determined by the cooling and escape processes of the secondaries, i.e., the 
pions, muons, kaons, and neutrons, whereas any additional cooling or escape processes of 
the primaries (protons, electrons, positrons) only affect the spectral shape and maximal 
energy. Therefore, flavor ratios may serve as the most robust prediction one may expect 
from an astrophysical source. However, they are much more difficult to measure: as mini- 
mal assumption, one needs to distinguish between muon tracks and cascades in a neutrino 
telescope. The measurement of cascades in this spirit has, for the first time, very recently 



been discussed in Ref. 14 for IceCube-22. 



Apart from the flavor ratio at the source, new physics in the neutrino propagation may 
also show up as another energy dependent effect. We have used two examples for such 
new physics effects, (invisible) neutrino decay and quantum decoherence, to illustrate the 
interplay between energy dependent flavor ratios and energy dependent new physics effects. 
Neutrino decay characteristically shows up at low energies, where the lifetime of the particles 
is shortest because of the smaller Lorentz boost. Quantum decoherence, however, may 
plausibly be a high energy effect, at least in certain scenarios discussed in the literature. 
While it is useful to have different source classes as a function of energy in a neutrino 
source to discriminate different decay scenarios, the considered quantum decoherence prefers 
sources other than a pure pion beam source. The reason is that both the standard case 
(flavor mixing only) and the quantum decoherence (even partial decoherence) case lead 
to approximate flavor equilibration for the pion beam case, which means that they are 
indistinguishable. 

As the final step, we have translated these requirements into our model parameters, where 
especially the magnetic field in the source plays an important role. For high magnetic fields, 
proton synchrotron losses limit the maximal energies, which can render high energy {n = 2) 
quantum decoherence invisible. In extreme cases, the neutrino energies may even be too low 
to pass the neutrino telescope threshold. For too low magnetic fields, the muon damping is 
not present in the flavor ratio which means that all decoherence scenarios and some decay 
scenarios are practically indistinguishable. Optimal magnetic fields are found to be in the 
range 10'^ Gauss ^ B < 10^ Gauss, where most effects are observable. Optimal (known) 
potential source classes may be AGN cores, white dwarfs, or GRBs, whereas AGN jets can 
only discriminate few effects as they typically perform as pion beams. 

We conclude that, although tough to measure, flavor ratios are a very interesting possibility 
to constrain particle physics properties using astrophysical sources in parameter ranges 
which would be otherwise not accessible. Instead of parameterizing the flavor composition 
at the source by constant (energy-independent) numbers, one could take into account the 
energy dependence of the flavor ratio which may be indirectly obtained from the multi- 
messenger connection. In many new physics scenarios, it is exactly the combination of this 
energy dependent flavor composition and the energy dependent new physics effect which may 
lead to new interesting observations or constraints - even with a single source, if enough 
statistics can be obtained. Although a high statistics measurement may not be expected 
from IceCube, future neutrino telescopes may be better optimized for cascade detection. 
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A Three flavor neutrino oscillation probability in presence of de- 
coherence 



Here we describe a phenomenological "bottom-up" framework to describe deviations from 
standard neutrino oscillation formalism in presence of decoherence irrespective of the details 
of the underlying theory {e.g., quantum gravity). Such a quantum system is described in 
terms of density matrices and the Liouville equation for the neutrino flavor density matrix 



p 54 - 56 , 74 - 78] (for any number of flavors) is modified to 

p = -i[H,p]+V[p] (26) 

where H is the Hamiltonian and p implies differentiation with respect to time. The first 
term is responsible for unitary evolution while the term T>[p\ contains effects due to non- 
unitary evolution. A commonly used form for T>[p] was given by Lindblad using quantum 
dynamical semi-groups [7^. For an A— level quantum system, it is possible to expand all 
the operators in the SU{N) Hermitian basis 

P = ^ [pM = \ [Po^ + Pi^i] 



1 
2 



H = ^ [h^,>^^,] = ^ [hoi + KXi 



V[p] = -[X^^d^uPu] (27) 



which leads to the equation of motion in the component form 

P^, = {h^u + d^y)p^ (28) 

where the subscripts /i, z/ depend upon the number of flavors, for three flavors, /x, z/ = 
0, 1, 8. The matrix elements h^^ are usually fixed from the form of the Hamiltonian while 
the elements of this matrix d^v can be fixed by assuming that the laws of thermodynamics 
hold. If one imposes the requirement of the monotonic increase of von-Neumann entropy 
(5* = — Tr[plnp]) which leads to hermiticity of operators, conservation of average value of 
energy, conservation of probability etc, we can further constrain the elements, d^y. For 
A = 3 case, the basis is given by eight Gell-Mann S'f/(3) matrices {Aj} along with the 



Identity /a matrix. We have the following matrix equation (from Eq. (28)) with eleven 
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decoherence parameters as in Ref. 54 , 
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where Uij 



1 



Smjj/4E with 6m'^j 
requiring Tr[p] = 1. We impose Tr[p] 



is the standard oscillation term. 



Po = 2/3 by 

1 and Tr[p] = which leads to conservation of 
probability. This leads to first row and column of the total matrix in above equation to be 
zero. Trace is preserved during evolution and hence po remains constant. The remaining 
eight components of p{t) can be obtained by solving the set of eight coupled equations 
subject to the initial condition 



PuM 



Ua2Ual Ua2Ua3 



(30) 



r/2 



and noting that the density matrix is given by 

I 3+P^+V3 
p{t) = - Pl+iP2 

\ P4 + ip5 



Pi - ip2 

-P3+1 
P6 + iP7 



Pa - iPb\ 
Pe - ip7 

2 2p8 , 

3 / 



(31) 



The diagonal elements of p{t) are referred to as populations while the off-diagonal elements 
as coherences. The phase information is contained in the coherences of the density matrix. 



Let us define the 8x8 block in Eq. (29) connecting the components pi{i 7^ 0) by £ 
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and then write the solutions pi{t) 


in 


terms of the elements of 



Pi(t) = pi(0)7Wn+p2(0)A^i2 + ...+P8(0)-M 
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P2(t) = PliO)M21+P2iO)M22 + ■.■+Ps{0)M28 



Psit) = Pi(0)7W8i+P2(0)A^82 + ...+P8(0)-M88 , (33) 



where Pi{0) are the components of the initial density matrix given in Eq. (30) and Eq. pT| . 
Finally the neutrino oscillation probability Pa^ can be computed using 

P^^it) = Tr[p,Jt)p,^(0)] , (34) 

where (0) is the "pure" neutrino density matrix corresponding to flavor Ujs at t = and 
Pu^t) is the density matrix at t for flavor i/q. 

For obtaining general expression of probability for astrophysical neutrinos, we need to av- 
erage over sin and cos terms, which gives us 

Po.p = \ + \iUl, - Ul2){Ul, - U%)D^ + \{Ul, + Ul, - 2Ul,){Ul + Ul, - 2U%)Ds , (35) 
where, Di^, = exp{— 2\l/t} and Dg = exp{—2St} are the decoherence-induced damping fac- 



tors. Thus out of the eleven decoherence parameters appearing in Eq. (29), only two (\l/ and 
6) appear in the flnal expression which are the A3 and Ag components in the decoherence 
matrix. If we look at the form of p{t), this implies that the coherences vanish for astrophys- 
ical neutrinos. This means that the rest of the physically allowed decoherence parameters 
are inaccessible by astrophysical neutrinos. If any of the phase information could be re- 
tained (which is what happens in atmospheric neutrino case) then the other decoherence 
parameters (corresponding to coherences in the p{t)) will also appear in the probability [54] . 
Let us now discuss the special cases in the large t limit: 



^ 7^ 0,5 = 0: 

1 1 
3 6 



Pap = :S + « (f^'i + - 2f/^3)(f^|i + U'p2 - 2f/|3) , (36) 



^ = 0,5 ^ 0: 



PaP = l + liUl, - Ul2){Ul, - UI2) , (37) 



^ = 6 = 0: 



PaP — \Uai\^ . (38) 

1=1 



$ = 5^0: 
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We note that only when both \1/ and 6 are non-zero, one gets the equihbrium steady state 

[so] . This means that a certain combination of non- 
zero parameters are responsible for obtaining equilibrium steady state answer or complete 
irreversibility. Violation of conservation of energy in the neutrino system is necessary but 
not a sufficient condition to ensure this. If one of the two parameters 6) are zero, the 
probability depends upon the mixing matrix elements in general. If both the parameters 
6) are zero, we recover standard oscillation result. We can refer to these three cases as 
achieving some sort of a steady state out of equilibrium. 



value for the probability (Eq. (39 
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